Fractional flux quanta at intrinsic metallic interfaces of 
noncentrosymmetric superconductors 
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We examine intrinsic interfaces separating crystalline twin domains of opposite spin-orbit coupling 
in a noncentrosymmetric superconductor such as CePtsSi. At these interfaces, low-energy Andreev 
bound states occur as a consequence of parity-mixed Cooper pairing, and a superconducting phase 
which violates time reversal symmetry can be realized. This provides an environment allowing flux 
lines with fractional flux quanta to be formed at the interface. Their presence could have strong 
implications on the flux creep behavior in such superconductors. 

PACS numbers: 74.20.Rp, 74.50. +r, 74.70.Tx 



Symmetry is a decisive factor for many properties of 
materials. Lowering a symmetry can yield new cou- 
plings between physical observables and causes intrigu- 
ing phenomena. The recently discovered noncentrosym- 
metric superconductors CePtaSi, CeRhSis, CelrSis, and 
Li2(Pt2;Pdi_a;)3B provide such examples [H, 0, BS]- 
these materials, the absence of an inversion center gen- 
erates antisymmetric spin-orbit interaction and leads, in 
the superconducting state, to parity-mixing of Cooper 
pairs, magnetoelectric effects, and many other interest- 
ing features li, 0, [§]■ many cases, such crystal 
structures permit the existence of twin domains exhibit- 
ing opposite inversion symmetry breaking within a single 
crystal. Actually, in the crystal growth processes of non- 
centrosymmetric materials, the formation of such twin 
domains is inevitable. The existence of twin domains 
in noncentrosymmetric superconductors is also suggested 
by a recent experiment, which revealed that a high qual- 
ity single crystal sample of CePtsSi exhibits a lower tran- 
sition temperature than polycrystal ones Since the 
origin of this behavior cannot be understood in terms of 
conventional impurity effects [lo| . possibly twin bound- 
aries could enhance the trend to superconductivity. Fur- 
thermore, recent NMR measurements of the single crys- 
tal sample are ingeniously interpreted by assuming the 
existence of twin domains (llj. Motivated by these ob- 



servations, in this letter, we investigate effects of intrin- 
sic interfaces between twin domains on the parity-mixed 
superconducting state. Our central finding is, that su- 
perconducting states with broken time-reversal symme- 
try can occur at the interfaces, allowing for fractional 
vortices. 

We consider a noncentrosymmetric superconductor 
such as CePtsSi and assume for simplicity a spheri- 
cal Fermi surface parametrized by the unit vector k = 
(cos Lp sin 9, sin if sin 0, cos 9) . The presence of a Rashba- 
type spin-orbit coupling, a(z x k) • s, induces a splitting of 
the electron bands and the Fermi surface into sheets, each 
exhibiting a specific spin structure. The superconducting 



phase displays a mixed parity [1,[1,[3|, and the state com- 
patible with experiments consists of an s- and a p-wave 
component, being of s ± p-character on the two Fermi 
sheets. Moreover, there is experimental evidence for a 
nodal gap structure, which suggests a dominant spin- 
triplet p-wave component with q — As/Ap < 1, where Ag 
and Ap denote the magnitudes of the s- and the p-wave 
components in the superconducting gap 12, T^]- For the 
calculations, we employ quasiclassical Eilenberger theory 
of superconductivity [H, 15, This method provides 
a convenient and powerful tool for describing supercon- 
ductivity and has been a ppl ied to noncentrosymmetric 
superconductors in Ref. |l3l|. According to Ref. [13], 
the superconducting state can be expressed by the so- 
called bulk coherence functions jb,Jb straightforwardly, 
which corresponds to the Riccati formulation of Eilen- 
berger theory 18, [l3|. Using an effective one-band de- 
scription, where the size of the band splitting is assumed 
to be small compared to the Fermi energy, we obtain [l3| 



IB = (Ty(7+(7+ +7_(T_), 

where the coefhcients are defined as 

A± 

7± = 
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(lb) 
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with A±(k) = As ± Apsm9 and ujn = (2n -I- VjirkBT 
denoting Matsubara frequencies. While we neglect the 
splitting of the bands, we keep their spin structure im- 
posed by the spin-orbit coupling, as described by 
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These spin matrices have the useful projection properties 
aj. — a± , (7+(7_ = (T_(T+ — and a+ + ~ 1. 

We now turn to the electronic properties of an inter- 
face separating regions A and B, which are characterized 
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B are then given by 



FIG. 1: Sketch of two intrinsic metallic interfaces in a noncen- 
trosymmetric superconductor with normal vector n perpen- 
dicular (left) or parallel (right) to the z-axis. a) The interface 
separates regions A and B, which exhibit a different direc- 
tion of the spin-orbit coupling a as indicated by the arrows, 
b) The situation is analogous to a junction between two or- 
dinary singlet superconductors with gap amplitudes A+ and 
— A_. The sign of the gap amplitude is illustrated by the two 
different colors. 



by the opposite sign of the antisymmetric spin orbit cou- 
phng according to a a — —as- An illustration of two 
specific situations for such an interface can be found in 
Fig. [U a), and similar setu ps h ave already been examined 

For our following analy- 



in a different context 
sis we neglect the direct influence of the interface on the 
superconducting order parameter and fix the moduli of 
the pair potentials to remain constant even at the inter- 
face. This approximation does not affect our discussion 
qualitatively, and quantitative corrections are minor. It 
leads to the simplification, however, that the coherence 
functions at the interface can be replaced by the corre- 
sponding bulk coherence functions. In the following, we 
use bulk coherence functions 7'^ and 7"* for region A ac- 
cording to Eqs. H]). Regarding region B, one might be 
tempted to get the bulk coherence functions 7^ and 7^ 
by simply interchanging the two gap amplitudes A±. For 
small values of q it is rather natural, however, to keep the 
dominant p-wave component Ap constant on both sides 
of the interface. Then, the s-wave component changes its 
sign across the interface and an additional phase factor 
— 1 has to be introduced on side B. We allow the gap 
function of region B to exhibit a further general phase 
difference with respect to region A, which is denoted by 
in the following. The bulk coherence functions of region 



1^ = {l-^++l+^-)(^ye"^ 
7^ = -(Tj,(7^CT+-f 7+'7-)e-**. 



(4a) 
(4b) 



Generally, once the coherence functions 7, 7 are known 
for a specific Fermi vector at a given point in space, also 
the quasiclassical Green's function g in 2 x 2 spin space 
is immediately available as 

5 = (1 - + 77) = 2(1 - 77)"' - 1- (5) 

The interface is implemented by well-established bound- 
ary conditions for the Green's function or the coherence 
functions, respectively [l^, [i^, [ll] . Restricting ourselves 
to a high-transparency interface, the resulting Green's 
function directly at the interface is given by 

g ^ 2{l~j^j^)-'-l 

= 2[1 - (7+<T+ + 7_<7_)(7_<7+ + 7+<7_)e-^^]-i - 1 

= 2[1 - (7+7-*^+ + 7-7+'^-)e"'1"' - 1 
2 

1, (6) 



1 - 7+7_e- 



where the projection properties of a± have been used in 
the intermediate steps. Note, that the expression Eq. 
([6]) for the quasiclassical Green's function at the inter- 
face only holds for quasiparticle trajectories with Fermi 
vectors k pointing from A to B. In the opposite case, the 
superscripts A and B have to be interchanged, and we 
find the symmetry relation ^(ciJn,— k) = g(ujn,i^)*- Sev- 
eral points should be mentioned here. Firstly, only bands 
having the same spin structure, a+ or contribute to 
this Green's function, the other combinations vanish by 
projection. Furthermore, the Green's function is propor- 
tional to the unit matrix. As a consequence, there is an 
analogy between this interface of noncentrosymmetric su- 
perconductors and a standard interface consisting of two 
singlet superconductors as illustrated in Fig. [TJ b). 

In the following, we derive the Josephson current den- 
sity through the interface. Using the symmetry relation 
stated above, we find 



j(0) = ineNoksTvp ^ (klm[g]). 



(7) 



a)„>0 



where (...) denotes averaging over half of the Fermi sphere 
determined by quasiparticle directions k pointing from 
region ^ to -B, and g is the unit matrix component of g 
according to Eq. Written in normalized quantities 

T = T/Tc and j = j/ATTeNoksTcVp, we eventually find 
the result 
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FIG. 2: (Color online) Part of the q, T-phase diagram of an 
intrinsic metallic junction in a noncentrosymmetric supercon- 
ductor. There are regions of a stable Josephson phase = 
(bottom) and (j> = tt (top). In between, a phase difference 
< (/) < TT is favorable. Here, the normal vector n of the 
interface is taken to be perpendicular to the z axis, corre- 



FIG. 3: (Color online) Part of the q, T-phase diagram of an 
intrinsic metallic junction in a noncentrosymmetric supercon- 
ductor. There are regions of a stable Josephson phase = 
(bottom) and = tt (top). In between, a phase difference 
< < TT is favorable. Here, the normal vector n of the in- 
terface is taken to be parallel to the z axis, corresponding to 



sponding to the lefthand scenario of Fig. [T] For T Tc, 0- the righthand scenario of Fig. [T] For T — > Tc, 0- and 7r-phase 



and TT-phase meet at q — V3/2 ~ 0.87. The dashed line in- 
dicates the boundary between 0- and vr-phase in the contrary 
scenario of a low-transparency tunnel junction. 



meet at q — l/\/2 ~ 0.71. The dashed line indicates the 
boundary between 0- and vr-phase in the contrary scenario of 
a low-transparency tunnel junction. 



where the upper (lower) formula stands for the situation 
with the normal vector n of the interface perpendicular 
(parallel) to the z-axis of the system. In both cases only 
components of the current flowing perpendicular to the 
interface are allowed by symmetry. The values of 7+ , 7_ 
in the integrand are real and depend on sin 9 themselves 
[cf. Eq. da)]. 

Numerical evaluation of the current-phase relations ac- 
cording to Eq. ([8]) allows us to determine the phase differ- 
ence (p of the stable interface states. We focus on > 0, 
keeping in mind that with (j) also —cj) is a stable solution. 
For the two situations depicted in Fig. [T] we can derive 
the phase diagrams displayed in Figs. [5] and [31 respec- 
tively. We find three regions in the g-T-phase diagram. 
In the region of small q, the stable state corresponds to 
(j) = and for q close to one it is = tt. The latter 
means, that the s-wave component would remain un- 
changed across the interface. Intriguing is a region in 
between these two limits, where the stable phase differ- 
ence has an intermediate value < < tt. Note, that the 
±0 solutions are degenerate for this intermediate region, 
reflecting the fact that such an interface state is time- 
reversal symmetry breaking, since changes sign under 
the time reversal operation. The intermediate region of 
q- values shrinks with increasing temperature, eventually 
reaching a single point at Tc. For the normal vector per- 
pendicular (parallel) to z this value is q = V3/2 ~ 0.87 
{q = I/V2 « 0.71) within our model. 

Qualitatively, the phase diagram can be understood 
as a consequence of Andreev bound states occuring at 
the interface. A sign-change of the gap function along 



a quasiparticle trajectory gives rise to the formation of 
zero-energy Andreev bound states. Such sign changes 
occur as soon as q > 0. Changing the phase difference 
from zero to a finite value may move these bound states 
away from zero energy, resulting in an energy gain ac- 
cordingly. For larger values of q the original spectral 
weight of the zero-energy Andreev bound states gets en- 
hanced. Consequently, upon increasing q a continous 
transition to a state of finite (j) occurs at some critical 
value qc{T). If q is increased further, eventually reaches 
the upper limiting value tt. 

The extent of the intermediate region in the phase dia- 
gram depends on the transparency of the intrinsic inter- 
face. The results presented sofar have been derived under 
the assumption of a high-transparency metallic junction. 
For comparison, we also examined the Josephson current 
for the opposite limit of a low-transparency tunnel junc- 
tion. Employing standard boundary conditions for the 

[illiil,!!!, we find 



quasiclassical propagators 



m=Dsin^-fY, / del -^j" 



where the following notation is used 

-1-1+ 



^ sin^ 9 c_L 

TT -I- 

sin 26 c\\ 



n, (9) 



=11 



arctan ^(lOa) 



(l + 7_7+)(7- -7+) l + 7_7+ 

-7-7+ 
(l-H7^)(l + 7i)' 



(10b) 



These results are valid to first order in the transparency 
D ^ 1, and, as in Eq. ([5]) for the metallic interface, the 
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upper formula corresponds to the orientation n _L z and 
the lower one to n || z. The main difference to the metal- 
lic case can be seen quite clearly: Since the current-phase 
relation in Eq. ^ is purely sinusoidal, the stable phase 
of the junction must be either or tt, depending on the 
sign of the amplitude factor. In particular, the interme- 
diate region has been shrunk to a single boundary line in 
the phase diagram. In Figs. [2]and[3l these boundaries 
between 0- and 7r-regions in the tunnel limit are sketched 
by the dashed lines for comparison. 

In the following, we concentrate on one remarkable 
physical consequence of the intermediate region where 

^ 0, TT. The degeneracy of the two phases ±0 gives rise 
to the possibility of line defects on the interface which 
carry fractional magnetic flux. They can exist at the 
interface only, and may generally exhibit fractional flux 
quanta $ according to 
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n G Z, 



(11) 



where $o — hc/2e is the standard flux quantum. As a 
consequence of this property, it is possible for a standard 
vortex to decay into two fractional ones on the interface, 
carrying the fractional flux (p/ir-^o and (1 — (j)/Tr) ■ $o, re- 
spectively. Both of these line defects are strongly pinned 
to the interface. If there are many of these fractional vor- 
tices lined up along the interface, they can act as a severe 
impediment for flux flow. Similar theoretical considera- 
tions have been made for domain walls in time reversal 
symmetry breaking superconductors 24 , 2^ 2^ . 

In summary, we find that interfaces between twin do- 
mains in a noncentrosymmctric superconductor such as 
CePtaSi could possess unusual properties. They can host 
low-energy Andreev bound states and, under certain con- 
ditions, give rise to a time reversal symmetry violating 
phase, a characteristic phase of the interface only. In this 
situation, fractional vortices could exist on the interface 
and severly influence the flux creep. Since the interface 
properties are different for different orientations, the flux 
creep properties would likely depend on the vortex direc- 
tion. Furthermore, the change of the phase cj) across such 
an interface can also modify special interference features 
of the Josephson effect in a magnetic field, if the inter- 
face intersects the junction between a noncentrosymmct- 
ric and a conventional superconductor. The low-energy 
Andreev bound states may be directly accessible by local 
tunneling probes such as scanning tunneling microscopes. 
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